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Abstract. We have searched for molecular gas towards 
the nucleus of four galaxies known to harbor a water vapor 
megamaser. CO(1->0) emission of NGC 2639 and NGC 
5506 was strong enough to allow us to map their inner 
regions. Weak emission from Mrk 1210 was detected and 
Mrk 1 was not detected at all. We report the tentative 
detection of the C0(2->1) line in NGC 5506. After this 
work, 12 of the 18 known galaxies harboring a water vapor 
megamaser have been observed in CO. 

The molecular gas content in the inner regions of water 
megamaser galaxies ranges from 5x 10 7 to 6x 10 9 Mq. The 
circumnuclear molecular gas surface density also extends 
over nearly two orders of magnitude. The maser luminos- 
ity is correlated neither with the total amount of molecular 
gas found in the inner few kpc of these galaxies nor with 
global properties of the molecular gas such as surface den- 
sity or filling factor; it is also independent of the infrared 
and optical luminosities. The only significant correlation 
we have found involves the maser luminosity and the low 
frequency radio continuum flux density. We conclude that 
the maser activity is intrinsically related to the energy of 
the active galactic nucleus whereas the intensity and even 
the presence of a water megamaser is independent of the 
molecular gas global properties such as the molecular gas 
content and surface density in the inner galactic regions. 

We have also found a possible anticorrelation between 
the molecular gas surface density and the rate of the mega- 
maser variations. A higher molecular gas abundance in the 
inner region could lead to higher maser variability because 
of larger nuclear flux variations due to the more variable 
gas infall, and/or because of more frequent interactions of 
the pumping agent with molecular gas condensations. 

Key words: galaxies: active - galaxies: Seyfert - galaxies: 
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1. Introduction 

The first water megamaser was discovered towards the 
nucleus of NGC 4945 (dos Santos & Lepine 1979). Strong 
emission at 22 GHz was detected, with a luminosity about 
100 times higher than that of W49, the most powerful 
galactic water maser. Until 1985 four additional water 
megamasers were discovered, towards the nuclei of the 
Circinus Galaxy, NGC 1068, NGC 4258 and NGC 3079. 
During the following decade no more extragalactic water 
megamasers were found. In 1994 a survey towards active 
galactic nuclei was carried out by Braatz et al. (1994), 
leading to the discovery of five new water megamasers, 
doubling the number known to that date. Recently, Braatz 
et al. (1996) have reported the discovery of 6 additional 
water megamasers towards active galactic nuclei. The last 
ones up to now have been found in NGC 5793 (Hagiwara 
et al. 1997) and in NGC 3735 (Greenhill et al. 1997b). 

Water megamasers, unlike the galactic and normal ex- 
tragalactic masers, have been detected at the nuclei of 
distant galaxies. Interferometric observations (Claussen & 
Lo 1986; Greenhill et al. 1995a; Greenhill et al. 1996) have 
shown that megamaser emission comes from within a re- 
gion around the galactic nucleus whose radius is in gen- 
eral smaller than 1 pc. Another important point is that all 
the galaxies harboring a water megamaser present some 
level of activity. This fact has led to suppose that the 
maser emission mechanism is identical to galactic masers. 
The huge difference in the energy involved is explained 
in terms of the pumping source: the energy source of the 
megamasers is the central object of the active nucleus, a 
much more powerful source than the central star powering 
the galactic masers. 

Braatz et al. (1997) have recently examined the condi- 
tions for detectability of water megamasers in terms of a 
variety of properties of the active galaxies, but not includ- 
ing the molecular gas content. It is known (Heckman et 
al. 1989) that Seyfert 2 have abnormally large quantities 
of dust and gas if compared with Seyfert 1 galaxies. On 
the other hand LINER galaxies are thought to be simply 
extensions of Seyfert 2 galaxies to lower luminosities, pho- 
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Table 1. Adopted parameters for the observed water megamaser galaxies 
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Obtained from the NASA Extragalactic Database (NED) 

6 See more information in the text (Sect. 2) 

c Obtained from the RC3 catalog (de Vaucouleurs et al. 1991) 



toionized by a weaker AGN spectrum (Osterbrock et al. 
1993). Thus, the fact that no water megamasers have been 
found in Seyfert 1 nuclei, while all of them are in either 
Scyfcrt 2 or LINER galaxies seems to indicate that a high 
concentration of molecular gas in the inner regions of the 
galaxies is a key parameter for the megamaser emission to 
be produced. 

To check this hypothesis and to find if the properties 
of the masers are related to those of the molecular gas, we 
have carried out a study of the molecular gas content and 
its properties in several of the water vapor megamasers 
for which no previous data existed or were incomplete. 
The megamasers we observed are four of those discovered 
by Braatz el al. in 1994: NGC 2639, NGC 5506, Mrk 1 
and Mrk 1210. In our analysis we have used the available 
CO data for other water megamaser galaxies found in the 
literature (Heckman et al. 1989; Planesas et al. 1989; Sahai 
et al. 1990; Aalto et al. 1991; Wang et al. 1992 and Young 
et al. 1995). 

2. Observations 

Most of the observations were carried out with the IRAM 
30 m radiotelescope at Pico Veleta (Spain) in July 1995. 
The transitions CO(l^O) at 115 GHz and CO(2^1) at 
230 GHz were observed simultaneously; the main beam 
diameters were 23" and 12" respectively, the main beam 
efficiencies, rjb, 0.68 and 0.40, and the forward efficiencies, 
r) fss , 0.92 and 0.86. The tentative CO(l^O) detection of 
Mrk 1210 obtained in this observing run prompted us to 
observe it again, observation that was carried out with 
the same instrument in May 1997, the antenna parameters 
being the same than in 1995. 

A 512-channcl filter-bank with a resolution of 1 MHz 
was used, providing a velocity coverage of 1200 km s -1 in 



the CO(1^0) line and 600 km s" 1 in the CO(2->l) line, 
wide enough to contain the whole range of the CO emis- 
sion. The velocity resolution was 1.3 km s -1 at 230 GHz 
and 2.6 km s -1 at 115 GHz. The focus was checked ev- 
ery few hours correcting it when necessary. Pointing was 
tested every 2 hours, only small errors were found, 2 — 6". 
The typical observing time for each individual point was 
2 hours. 

All observations were made nutating the subreflector 
at a rate of 1 Hz, to a position 4' in azimuth from the ob- 
served point. Thus, we were able to obtain the flat base- 
lines required for the detection of the wide lines expected 
towards the inner regions of these galaxies. Linear base- 
lines have been substracted from each spectrum. To obtain 
the absolute calibration of the temperature scale, spectra 
towards IRC+10216 and Orion IRC2 at the redshifted fre- 
quencies were taken and compared with the IRAM 30 m 
catalog of calibrated molecular lines (Mauersberger et al. 
1989). The temperature scale in figures is the antenna tem- 
perature corrected for atmospheric attenuation and rear 
spillover, Tj£. The CO luminosities and derived parame- 
ters have been converted to the TJj temperature to allow 
compatibility with the data obtained from the literature. 

The adopted parameters for the galaxies we observed 
are shown in table 1. Distances were determined from 
Vlsr and a Hubble constant of 75 km s _1 Mpc -1 . Infrared 
and far-infrared luminosities were computed using the ex- 
pressions given in Perault (1987). Isotropic maser lumi- 
nosities are given in table 2. Due to the strong tempo- 
ral variability of the emission, they are only valid for the 
epoch its observation was carried out. 

3. Observational results 
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Fig. 1. CO(l^O) emission profiles towards the central region (23" diameter) of the three detected water megamaser galaxies. 
Velocity resolution is 21 kms -1 . The two NGC objects show double peak spectra, which may indicate the existence of a molecular 
gas ring surrounding the central AGN at a kpc scale. The spectrum of Mrk 1210 is single-peaked, a possible indication that the 
expected molecular gas ring is seen face-on. This spectrum was obtained in May 1997 



3.1. CO (1^>0) line 

The spectra obtained towards the central regions of the 
three detected galaxies are shown in Fig. 1. No emission 
was detected towards Mrk 1 down to a rms noise of 10 
mK for a resolution of 42 kms -1 . The spectra of NGC 
2639 and NGC 5506 are both double peaked. In NGC 2639 
the centroid of the CO emission is at Vlsr = 3258 kms -1 
rather similar to the LSR velocities derived from optical 
observations (3163 kms -1 ) or from the HI 21 cm line ob- 
servations (3298 kms -1 ). In NGC 5506 the LSR velocity is 
1827 kms -1 from optical observations, 1813 kms -1 from 
HI 21 cm line observations whereas the CO (1-^0) centroid 
velocity is about 1880 kms -1 . In Mrk 1210 the centroid 
velocity of the single fitted Gaussian function appears at 
Vlsr — 3995 kms -1 in good agreement with values deter- 
mined from optical observations (3968 kms -1 ) or from the 
HI 21 cm line observations (3978 kms -1 ). The intensity is 
20 mK and the rms noise is 4 mK for a resolution of 21 
kms -1 . 

Small maps of NGC 2639 and NGC 5506 were obtained 
(Fig. 2). The typical rms noise for the thirteen NGC 2639 
spectra is between 2 and 7 mK. The mapped region has a 
diameter of 40" which corresponds to a size of 8 kpc. The 
emission is extended in the direction NW-SE, in agree- 
ment with the position angle of the galaxy (140°). The 
distribution of the integrated emission has an elliptical 
shape with major and minor axis of 30x16" or 6x3.2 kpc, 
in agreement with the 58° inclination angle of the galaxy. 

For NGC 5506, in addition to the center we observed 
four positions 20" away from the nucleus. The rms noise 
lies between 4 and 6 mK. The observed region is about 
60" in diameter which corresponds to a size of about 7 
kpc at the source distance. The emission is more extended 
towards the E-W direction, in agreement with the position 
angle of the galaxy (91°). 



3.2. CO (2-^1) line 

Observations of the CO(2— >1) line were carried out simul- 
taneously to those of the CO(1^0) line. We obtained a 
13-point map of NGC 2639, a 5-point map of NGC 5506 
and single central spectra in Mrk 1 and Mrk 1210. The 
rms noise achieved was 4 mK in the NGC 2639 and Mrk 
1210 spectra and 40 mK in that of Mrk 1. 

We detected tentatively only the central region of NGC 
5506. The rms noise, for a resolution of 21 kms -1 , is 8 
mK and the line integrated intensity 5.5 K kms -1 . The 
centroid of the CO(2 — >1) emission is at 1870 kms -1 . In 
Fig. 3 (right panel) the CO (2— >1) spectrum produced by 
adding the five spectra obtained for this galaxy is shown. 
The CO(2-»1)/CO(1-»0) line ratio in the inner 20" region 
is 0.72, showing that the gas is optically thick in 12 CO and 
that the excitation temperature is about 8 K (cf. Braine 
& Combes, 1992). 

3.3. CO luminosities 

The CO luminosities in table 3 have been obtained as- 
suming a Gaussian antenna beam, using the expression: 
L co = 1.13 d\ I C o, where J CO (K(T^)kms -1 ) is the in- 
tegrated intensity derived from the CO(l— >0) spectra and 
db (pc) is the antenna beam width at the source distance. 

The H2 mass has been derived from the CO luminos- 
ity using the standard conversion factor found in galactic 
giant molecular clouds, Mh 2 = 5.8 Lco- High gas tem- 
perature and high metalicity could lead to an overestimate 
while high density and photodissociation of CO could lead 
to an underestimate of the amount of molecular gas, so, we 
have used the standard conversion as the most adequate 
choice. 

The molecular gas surface density, Sh 2 has been deter- 
mined for the region of the galaxy included in the antenna 
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Fig. 2. CO(l— »0) maps of the circumnuclear regions of NGC 2639 and NGC 5506. In NGC 2639 the velocity range extends 
from 2700 to 3700 km s -1 and in NGC 5506 from 1450 to 2300 kms -1 . The resolution of the spectra is 21kms~ 1 in NGC 5506 
and 42 kms -1 in NGC 2639 



beam, except in the galaxies where this size is larger (see 
table 2) than the CO effective diameter (Z?co otf > defined as 
the diameter that contains the 70 % of the total CO emis- 
sion). In the last case we used the CO effective diameter 
to calculate the surface density of molecular gas. To de- 
termine this size we used the optical diameter D25. From 
the data of Kenney & Young (1988) where CO effective 
diameters are given for a set of galaxies, we have found a 
correlation between -Dco eff and -D25 that is roughly inde- 
pendent of the model chosen for the CO emission distri- 
bution. The relation is: D COoff = 0.44 D25, very similar to 
that obtained by Young & Scoville (1991). 

In table 2 the parameters related to the molecular gas 
are shown for the four galaxies we observed as well as 
for other water megamaser galaxies whose data have been 
taken from the literature. 

4. Results 

4.1. Mrk 1210 

Our data show that Mrk 1210 is a peculiar object among 
water megamaser galaxies. It has a relatively low surface 
density of molecular gas, not much higher than that of 
normal galaxies like the Milky Way or M33. However, its 
star formation efficiency, estimated from the Lir/Mh 2 ra- 
tio, is the highest in our sample. The reason may be that 
in this kind of objects, a significant fraction (up to 50 %) 
of the IR luminosity could come from the active nucleus 
itself, as is the case for NGC 1068, resulting in an over- 



estimate of the star formation efficiency up to a factor of 
two. 

Nevertheless, there are other indications of a high 
level of recent star formation activity in this galaxy. The 
S25/S100 ratio (where S25 and sioo are the IRAS fluxes) has 
been suggested as the best tracer of the ionizing photon's 
source (Dultzin-Hacyan ct al. 1990). The high measured 
ratio, S25/sioo=l-59, in this galaxy would thus indicate a 
large star formation activity, in spite of the shortage of 
molecular gas. 

4-2. Molecular gas properties and megamaser variability 

One of the most interesting properties of water mega- 
masers is their variability. The observations carried out up 
to now have shown that the megamaser features vary in 
strength as well as in velocity (e. g. Claussen & Lo, 1986; 
Greenhill et al. 1995b; Baan & Haschick 1996; Braatz et 
al. 1996; Hagiwara et al. 1997). These intensity variations 
can reach a factor of ten. The time scale of the variations 
range between some weeks and a few years, though Green- 
hill et al. (1997a) have recently found intensity variations 
in the Circinus galaxy on a time scale of a few minutes. 

We have tried to determine whether long-term time 
variations could be related somehow to the abundance of 
molecular gas. The standard AGN model assumes that 
Seyfert 2 galaxies have a parsec-scale thick torus of gas and 
dust that blocks the direct view of the nucleus (Antonucci 
& Miller 1985). In addition to that, a large amount of 
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Fig. 3. Left and center panels: CO(l— >0) integrated emission of the circumnuclear regions of NGC 2639 and NGC 5506. The 
spectra have been obtained by adding the individual spectra of all the observed points (5 in the case of NGC 5506 and 13 
in NGC 2639). Right panel: CO(2— »1) integrated emission of the circumnuclear region of NGC 5506, it has been obtained by 
adding the spectra of the five observed points. The spectral resolution is 21 kms -1 in all cases. The intensity scale range is also 
the same, to allow the visual comparison of the intensities. The velocity dispersion in NGC 2639 is very high, the CO(l— >0) 
line has a width of 600 kms" 1 . In NGC 5506, the emission appears whitin the same range of velocities (1700-2100 kms -1 ) for 
both observed transitions 



dense molecular gas has been found in the inner regions 
of these galaxies (cf. Planesas et al. 1997). 

The distribution of this gas could be relevant to the 
diverse phenomena produced by the interaction of the ra- 
diation or matter ejected from the nuclear source. The 
megamaser phenomenon could arise in the external lay- 
ers of the nuclear torus, its long-term variability linked to 
the inhomogencities of the gas distribution in the exter- 
nal layers of the torus, but could also be affected by other 
structures that are known to exist in the inner galaxy. In 
fact, a study of the circumnuclear structure of the galaxy 
NGC 3079 (Baan & Irwin 1995) showed that there are 
several structures of molecular gas at different radii from 
the nucleus. 

On the one hand the observed megamaser variations 
could be related to the intrinsic nuclear variations. Intrin- 
sic variations at different time scales in radio-continuum 
flux density have been found in one Seyfert 2 galaxy 
(NGC 1275, Nesterov et al 1995): slow variations have 
typical time scales of 20-30 years, faster variations hav- 
ing time scales ranging from some months at millimeter 
wavelengths to a few years at lower frequencies. All the 
radio-continuum variations are thought to be produced 
by changes in the accretion rate of material towards the 
central black hole. In this case the molecular gas avail- 
able near the nuclear region (scale of kpc) fuels the ac- 
cretion disk of the AGN and any variation in the falling 
rate of gas towards the disk could eventually produce the 
observed variation in the nuclear radio-continuum emis- 
sion. The gas present at kpc scales can be transferred to 
the nuclear region by a kpc-size bar that seems to be very 
common in spirals (including Seyferts), as shown by NIR 



observations (Mulchaey et al. 1997). Hence, the presence 
of molecular gas at a scale of hundreds of parsecs from 
the nucleus could be significant in the generation of the 
observed variability of the radio-continuum and therefore 
of the maser intensity. 

Another contribution to the observed megamaser vari- 
ability could arise from the small structure in the parsec- 
sized torus or in Giant Molecular Clouds (GMCs) located 
not far from the nuclear source. In particular, when one 
of such clouds crosses our line of sight towards the central 
maser emitting region it could act as an extra amplifier of 
the emission. The maser variability might be due in part 
to the inhomogencity of the molecular clouds, those re- 
gions with higher density would produce an increment in 
the observed maser intensity. Taking a maximum size of 
800 AU (Falgarone et al. 1992) for the gas accumulations 
inside a cloud, and assuming a rotation curve similar to 
that of NGC 891 (e.g. Garcia-Burillo et al. 1992), the time 
scale expected for variability produced in this way would 
be less than a few tens of years, compatible with the long 
period variability observed in megamasers. Accepting this 
reasoning, the 3 year average time-scale of the variability 
found in water megamasers (cf. Table 2) would be pro- 
duced by structures ten times smaller. 

It could be expected that in both high abun- 

dance of molecular gas in the inner region could lead 
to a higher variability of the maser, via variable infall 
rate to the nucleus or via the interaction of the pump- 
ing source with inhomogeneous ambient clouds. One of 
our results agrees with this hypothesis: we have found 
an anticorrelation (Fig. 4, center panel) between Sh 2 m 
the central region of water megamaser galaxies and the 
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Table 2. Molecular gas and maser parameters in the inner region of the water megamaser galaxies 
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7.4 3 
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NGC 4945 
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19.6 
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9.7 


270 
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57 
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8.2 4 
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3.13 


11.7 


24 


NGC 1386 


2.3 
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7.2 5 


8.0 


30 
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120 


NGC 5347 


8.9 
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8.3 6 
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2.25 


32 
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7.6 


7.0 7 
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38 


16 
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a Linear size of the antenna beam at the galaxy distance 

6 Effective CO diameter, defined as the diameter that contains 70% of the CO emission (see more details in Sect. 3.3) 
c CO luminosity in the central region of diameter 0b- Obtained from the CO(l— >0) line. Reference numbers correspond to: (1) 
This paper; (2) Planesas et al. 1989; (3) Young et al. 1995; (4) Aalto et al. 1991; (5) Sahai et al. 1990; (6) Heckman et al. 1989; 
(7) Wang et al. 1992 

d Obtained using the standard conversion factor found in galactic GMCs, Mh 2 = 5.8 Leo 

e Surface density of molecular gas. It has been obtained for a region whose diameter is the smallest value between 9 b , and 
^cOeff (see details in Sect. 3.3) 

1 Isotropic maser luminosities, taken from Braatz et al. (1996) 



rate of the relative variations of the megamaser intensity: 
(-fmax/-7min)/At, where 7 max (7 mi „) is the maximum (min- 
imum) peak flux density observed in the megamaser at 
22 GHz and At is the time elapsed between them. The 
main objection that can be raised to this tentative result 
is the low spatial resolution of our determination of £h 2 , 
as our measurements include a more extended region than 
the expected one relevant for the variability of the nuclear 
source. 

4-3. Molecular gas properties and megamaser intensity 

From our estimates of £h 2 we have calculated the number 
of clouds present in the region we observed, using as the 
typical GMC parameters a diameter of 40 pc and a mass 
of 3xl0 5 M Q (cf. Wilson et al. 1990). The total volume 
of the clouds gives us the filling factor in the observed 
region. Our results show that there is no correlation be- 
tween the isotropic L mas0 r and the filling factor, see Fig. 4 
(left panel) . This implies that the total luminosity of the 
megamaser is not affected by the abundance of molecu- 
lar clouds in the inner kpc, what is reasonable taking into 
account that the typical filling factor is ~ 0.01. 

Another naively expected correlation that is not found 
from our data is that of log(L masor ) vs. S H2 - If we suppose 
that the optical depth of a molecular region is propor- 
tional to Sh 2 and that the amplification is exponential: 
I = I cxpr (Baan & Haschick 1996), a dependence be- 
tween Eh 2 and log(L masor ) is expected. However, as can be 



seen in the Fig. 4 (right panel) there is no correlation. Due 
to the complex medium the radiation has to go through 
it would not be strange that global amplification is far 
from exponential. Moreover, most of the amplification is 
likely to take place in a much smaller scale than the one 
considered when estimating Sh 2 • 

4-4- Energy of the central source and maser luminosity 

We have not found any correlation (see Fig. 5) between 
L maser and any other relevant luminosity (i.e., £fir.) -^IR) 
Lb, Lx)- This is surprising, as a closer relationship be- 
tween .Lfir and L x with L masor could be expected. On the 
one hand, the FIR radiation is thought to be re-emission 
by molecular and dust clouds of the high energy radiation 
coming from the central source. On the other hand, the 
adequate pumping conditions are thought to be produced 
by the heating of the clouds by X-ray shocks. Therefore, 
although one would expect the maser luminosity to be re- 
lated to the energy release by the central source, this is 
not the case. 

A better measure of the intensity of the central source 
could be the radio continuum emission. We searched for a 
correlation between the maser luminosity and the absolute 
low frequency radio continuum flux density of the galaxies. 
The radio continuum emission is the sum of nuclear source 
emission plus the extended emission associated with star 
formation; however, most of the radio continuum emission 
in active galaxies comes from the nucleus. Braatz et al. 
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Fig. 4. Maser characteristics vs. molecular gas properties. Left panel: L mascl vs. molecular gas filling factor of GMCs (see 
Sect. 4.3). To estimate the volume filling factor we have considered two cases: (i) a cylindrical geometry (represented in the 
figure by the black circles) with a diameter the size of the antenna beam at the source distance and a thickness of 300 pc (cf. 
Garcfa-Burillo et al. f 992) and (ii) a spherical geometry (represented by the empty circles) with the diameter of the antenna 
beam at the source distance. The filling factor is a measure of the probability that a molecular cloud is in the line of sight of 
the maser. No correlation is found, therefore the isotropic luminosity of the megamaser seems to be independent of the amount 
of molecular gas present in the central region of the galaxy. L mascr is given in Lq. Center panel: (7max//min)/A£ vs. £h 2 ; the 
first parameter is a measure of the variation rate of the megamaser intensity. There is a possible correlation between the two 
magnitudes (correlation coefficient of -0.71). Time is given in months and Sh 2 in M pc~ 2 . This result may be affected by the 
fact that the monitoring of the last discovered megamasers is still incomplete. Right panel: log(L mase r) (in L Q ) vs. Eh 2 (in 
M Q pc -2 ). No correlation is found, see Sect. 4.3 



(1997) have found a dependence of the detectability of 
water megamasers on intrinsic nuclear radio intensity at 
cm wavelengths. In the lower panels of Fig. 5 £ masor is 
plotted vs. the radio continuum flux data available for the 
known water megamaser galaxies. Weak correlations are 
found between £ masor and the flux density at 1.4 and 8.4 
GHz, although the number of data points at 8.4 GHz is 
quite small. 

From the correlations between L masor and the radio 
continuum emission at 1.4 and 8.4 GHz may be inferred 
that the maser luminosity is an intrinsic, energy-related 
property of the galactic nucleus (which could be charac- 
terized by the radio luminosity) rather than by the inner 
galactic regions (probably best characterized, e.g., by the 
infrared luminosities and molecular gas abundances). 

5. Discussion 

We have found no correlation between the maser lumi- 
nosity and other luminosities, from far-infrared to X-ray 
luminosity, for a sample of a dozen galaxies, in spite of 
the presumption that those luminosities could be related 
in one way or another to the generation or pumping of the 
maser (cf. Sect. 4.4 and Fig. 5). This result suggests that 
the maser phenomenon is very localized and is not related 
to global properties (mass or luminosity) of the galaxies. 
This is further supported by the correlation found between 
the low frequency radio continuum and the maser luminos- 
ity, both arising mostly very close to the galactic nucleus. 

The accumulation of clouds of dense molecular gas 
around the nucleus of these galaxies seems necessary for 



the water megamaser to be produced. However, we have 
found no apparent relation between the maser luminosity 
and the surface density or the total content of molecular 
gas in the inner galactic regions (say, the innermost kpc). 
In spite of that, molecular gas is known to strongly shape 
the appearance of the central regions of galaxies at dif- 
ferent wavelengths. High angular resolution (Martin et al. 
1989) and high sensitivity CO maps (Krause et al. 1990) 
of the water megamaser galaxy NGC 4258 have shown the 
importance of the molecular gas distribution in confining 
the radio jets arising from the nucleus. Cecil et al. (1995) 
have shown the coincidence of the radio jets and soft X- 
ray jets in this galaxy; the hot gas producing the thermal 
X-ray jets was interpreted by these authors as shocked gas 
pulled away from the ambient molecular clouds. 

The only possible correlation we have found (Fig. 4, 
center panel) involves the molecular gas surface density in 
the inner regions of the galaxies, in the sense that it is an- 
ticorrelated with the variation rate of the maser intensity. 
This result can be interpreted in terms of the maser emis- 
sion being produced or further enhanced by the interaction 
of the nuclear jets with clouds of matter surrounding the 
active nucleus at different scales. The combined effect of 
cloud movement, the small scale structure of the clouds 
and the intrinsic variability of the central source (itself re- 
lated to the gas infall towards the nucleus) could account 
for the observed variability of the megamaser intensity. 
This scenario is neither meant to explain maser variability 
at all time scales, nor all the subtleties observed in water 
maser variability, but to qualitatively show that one as- 
pect of the variability (its rate of amplitude variations) 
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Fig. 5. Plot of -Lmascr vs. several relevant parameters for the water megamaser galaxies. The units of Lir, Lfir and L maser arc 
Lq, H2 mass is in M . Radio fluxes are given in mjy, X-ray fluxes in erg cm -2 s _1 and distances in Mpc. There is a certain 
correlation between the isotropic maser luminosity and the radio continuum flux density at 1.4 GHz and 8.4 GHz (correlation 
coefficients of 0.6 and 0.87 respectively), see Sect. 4.4 



could be somehow related to the abundance of molecular 
gas around the nucleus. 

6. Conclusions 

We have searched for molecular gas towards the nucleus of 
four galaxies known to harbor a water vapor megamaser, 
and detected the CO(l— >0) emission in three of them and 
the CO(2— >1) emission in one. With this work 12 of the 
18 known water megamaser galaxies have been observed 
in CO, and only the most distant of the observed ones, 
Mrk 1, has not been detected yet. 

1. The 12 water megamaser galaxies with molecular gas 
data available arc not an homogeneous set regard- 
ing their molecular gas properties. The amount of H2 
in their circumnuclear regions ranges from 5xl0 7 to 
6xl0 9 M Q . The extreme values of the H2 surface den- 
sity, Sh 2 , for the central kpc are 6 and 280 M pc~ 2 . 
This parameter extends over a range of 2 orders of 
magnitude, a range similar to that of Seyfert galaxies, 
starburst galaxies or luminous infrared galaxies. The 
maser luminosity, i masor , is not correlated to the to- 
tal molecular gas mass. Therefore it seems that the 
total amount of molecular gas in the inner few kpc is 
not a fundamental parameter on which depends the 
existence and the average intensity of the water mega- 
maser. 

2. L mascr is not correlated with Eh 2 or to the filling fac- 
tor of giant molecular clouds. Apparently the maser 



luminosity does not depend on the content of molecu- 
lar gas in the inner kpc. However, the accumulation of 
clouds of dense molecular gas is believed to be neces- 
sary for the generation of a water megamaser. Obser- 
vations with higher angular resolution of the molecular 
gas in the inner regions would help to solve the issue. 

3. The only correlation we have found involving the maser 
emission and molecular gas parameters is between the 
rate of relative variation of the maser intensity and 
Eh 2 • This fact may indicate that a high abundance of 
molecular gas in the inner regions could lead to higher 
variability in the maser emission, on the one hand, due 
to the higher variability of the central pumping source 
produced by wider variations in the gas infall; on the 
other hand, due to the more frequent interactions of 
the pumping agent with molecular gas condensations. 

4. imasor is not correlated to any other luminosity (in- 
frared, optical, X-ray, blue). However, we have found 
some correlation between L mascr and the global ra- 
dio continuum flux density at 1.4 and 8.4 GHz. This 
fact supports the idea that L maS cr is a property re- 
lated to the galactic nucleus (characterized by the ra- 
dio luminosities) rather than to the inner galactic re- 
gions (characterized by the infrared luminosity and the 
molecular gas content). 

5. Mrk 1210 stands as a peculiar object, having the high- 
est star formation efficiency among water megamaser 
galaxies is spite of its relatively low molecular gas con- 
tent. 
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